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ABSTRACT 

We present a detailed analysis of the spatial distribution of the stellar populations 
in the Galactic globular c luster ui Centauri. Ta king advantage of the large photomet- 
ric catalog published bv iPancino et al.1 l)2000|h we confirm that metal- rich popula- 
tions have a spatial distribution which is significantly different from the metal-poor 
dominant population. In particular: (i) the different sub-populations have different 
centroids and (ii) the metal-poor population is elongated along the E-W direction, 
while the metal-rich populations are oriented along the orthogonal direction, i.e., N-S. 
The evidence pr esented here ca n partially explain the weird spatial metallicity segre- 
gation found bv lJurcsIkl 1^998), and further supports the hyphothesis that different 
sub-populations in oj Centauri might have had different origins. 

Key words: globular clusters: individual: NGC 5139 



1 INTRODUCTION 

The peculiar nature of the Galactic globular cluster ui Cen- 
tauri has been known and studied for more than forty years. 
Besides being the most massive and luminous in the Milky 
Way, it is presently the only globular cluster that shows a 
spread in its heavy elements content. Recent findings have 
shown that (i) at least three primary enrichment peaks do 
exist in this clust er, including the recently discovered m etal- 
rich component jLee et al.l ll99Sl : IPancino et alJlioOOtl : (ii) 
an age spread of 3-5 Gyrs see ms to be required to explain 
the turnoff region morp hology jHughes fe WallersteinllioOOt 
Hilke r^g Richtlerl2000t) : (Hi) the chemical enrichment of the 
metal-poor and intermediate stars is mainly due to the re- 
tention of SNe II and intermediate mass A GB stars ejecta 
JSmith et al.l200d : lNorris fc Da Costall995lh (iv) the metal- 
rich stars appear to have a low er a-enhancement, m ost prob- 
ably due to SNe la pollution iPancino et al.ll2002lh 

All these pieces of evidence suggest that u Centauri 
could be an intermediate object between normal globu- 
lar clusters, which are unable to retain any of the super- 
novae ejecta, and the dwarf spehroidal galaxies (dSph), 
which are the smallest stellar systems capable of self- 
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enrichment, u Centauri could also be the remaining nucleus 
of a dwarf galaxy that was stripped during its interaction 
with the Milky Way, in possible analogy with the complex 
and still debated case of M 54 and the Sagittarius dSph. 
The possibility that uo Centauri comes from "outside" the 
Milky Way seems al so required to explain its present orbit 
jPinescu et al]ll999lh 

However, other clues complicate the picture. For ex- 
ample, the unusually high ellipticity of u Cen, that 
has been demonstrated to be sustained by rotation 
jMerritt fc Tremblavl |l994h . is compatible with the flat- 
tene d shapes resulting from the merger of two globular clus- 
ters iMakino. Akivam a fc Sugimqto]| l99llh and anyway the 
long relaxation time lDk^rgoTOkil^993! : iMerritt et al.lll997l) 
grants that u> Centauri is not completely relaxed dynami- 
cally. Moreover, Norris et al. (1997) showed that only stars 
with [Fe/H]^ —1.2 in uo Cen do rotate, while the more metal- 
rich components show no evident sign of rotation. Pancino 
et al. (2000) showed that, while the metal-poor population 
exhibits the well known E-W elongation, the two metal-rich 
populations show a more pronounced ellipticity, but with an 
elongation in the N-S direction. Finally, Ferraro et al. (2002) 
have shown how the most metal-rich population shows a dif- 
ferent bulk proper motion with respect to the whole cluster. 
These pieces of evidence point toward a different dynami- 
cal origin for the various sub-populations in uj Cen, possibly 
resulting from a major merger event in the cluster's past 
history. 
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Figure 1. Isodensity contour lines for the three RGB samples 
denned in the text: the RGB-MP (upper panel), the RGB-Mint 
(middle panel) and the RGB-a (lower panel). The contour lines 
are normalized to their peak density and are plotted in steps of 
10% of the peak density. In all panels the axes show the distance 
from the cluster centre (700±20 pix, 1900±20 pix), in arcminutcs. 



Given this framework, we have started a coordi- 
nated spectro-photometric study of to Centauri (see 
Pancino et all l200fl 120021: iBellazzini. Ferraro fc Pancinol 
20011: iFerraro. Pancino fc Bellazzinil 2002aF 

Fterraro! Bellazzm^^Pancmon2002bilc)riglia et~al I 12003ft . 
specifically devoted to the characterization of the sub- 
populations, and aimed at understanding the origin and 
evolution of this complex stellar system. In this paper, 
we exploit the large number of stars (more than 230,000 
in total) and the wide area coverage (33' x 33') of the 
photometry published bv lPancino et all (l2000l) . to study in 
deeper detail the structural properties of the red giants in 
lu Centauri. 

The article is organized as follows. In Section [5] we 
define the photometric samples corresponding to the sub- 
populations of giants in uo Centauri; in Section [3] we study 
the surface density distributions and concentration of the 
three sub-samples; in Section 2] we measure the ellipticity, 
the centres positions and the orientation of the three sub- 
populations; in ScctionElw e comm ent on the metallicity seg- 
regation found bv I Jurcsikl il998T> . Finally, in Section |SJ we 
summarize our main results and discuss them in the frame- 
work of the present theories about the formation and evolu- 
tion of uj Centauri. 



2 SAMPLES DEFINITION 

Using the metall icity informat i on fro m the low resolution Ca 
triplet survey bv lNorris et alJ (Il99d) and the morphology of 
the RGB from the colour-magnitude diagram, Pancino et 
al. (2000, see their Fig. 2) identified three sub-populations 
of RGB stars, with different average metallicity and photo- 
metric properties. To analyse in detail the structural prop- 
erties of these sub-p opulations, we have extracted from the 
IPancino et alJ ifcOOOT) catalogue the following photometric 
samples (with B < 16 mag): 

(i) the RGB-MP sample, corresponding to the main, 
metal-poor peak of the metallicity distribution, around 
[Ca/H]~-1.4 or [Fe/H]~-1.7. This population comprises 
~ 70% of the whole RGB population, and our photomet- 
ric sample contains 2630 stars. 

(ii) the RGB-Mint sample, comprising the secondary, 
intermediate metallicity peak around [Ca/H]~-1.0 or 
[Fe/H]~-1.2, together with the long, extended tail to higher 
metallicities. This sub-population accounts for ~ 25% of the 
whole RGB population, and our sample contains 816 stars. 

(iii) the RGB-a sample, the newly discovered metal-rich 
population that comprises ~ 5% of the RGB stars and has a 
meta llicity of [Ca/H]~-0.5 and [Fe/H]~-0.6 JPancino et alJ 
2002). Our photometric sample contains 128 stars. Although 
this last sub-population has many fewer stars, it is the most 
numerous sample presently at our disposal, and its size could 
be significantly increased only when the RGB-a counterparts 
in other evolutionary phases (i.e., horizontal branch, sub- 
giant brach and main sequence) will be identified. 



3 SURFACE DENSITY DISTRIBUTIONS 

We computed the surface density distributions of the three 
sub-populations defined in Section [5] using a fixed kernel 
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Table 1. Results of the two-dimensional generalization of the 
K-S test for the three sub-populations. The first column shows 
the populations that are actually compared, the second shows 
the maximum difference in the cumulative distributions D, while 
the third column shows the derived probability P that the two 
populations are drawn from the same parent distribution. 

Population D P 

RGB-a vs. RGB-MP 0.171 0.012 

RGB-Mint vs. RGB-MP 0.068 0.039 
RGB-Mint vs. RGB-a 0.117 0.225 



Table 2. For each sub-population (column 1), the equivalent ra- 
dius, in pixels, of the ellipses used to approximate the 90% iso- 
density contour lines (close to the peak) is shown in column 2 
( r e 90%) an d its ratio with the radius of the RGB-MP population 
in column 3 (Rgo%)- Columns 4 and 5 show the corresponding 
values (r e qq% and Rgo%) f° r the ellipses used to approximate the 
60% isodensity contour level (close to half maximum). 

Population r e90% R 90% r e60% R 60% 

RGB-MP 213.8±6.5 1.00 578.9±13.5 1.00 
RGB-Mint 169.3±15.9 0.79 601.7±7.6 1.04 
RGB-a 66.2±5.5 0.31 323.1±5.0 0.56 



estimator algorithm dSilverman||l986t iMerritt fc Tremblavl 
ll994l:ISeleznevll99Sl) . with a kernel half- width of 500 pixels 1 
and a grid of 100 pixel cells. The resulting surface density 
plots are shown in Figure where the isodensity contour 
lines shown are normalized to the maximum density of each 
distribution, in steps of 10%. From now on, we will refer to 
each isodensity level using the fraction of the peak value, i.e., 
0.3 for the isodensity level corresponding to 30% of the peak 
value. The peak density values are: 48.3 stars arcmin -2 for 
the RGB-MP, 14.6 stars arcmin~ 2 for the RGB-Mint and 
2.3 stars arcmin -2 for the RGB-a. 

A first, qualitative comparison of the three distributions 
shown in Figure reveals the following general facts: 

(i) The RGB-MP population is clearly elongated along 
the E-W direction, reflecting the well known elliptical shape 
of the whole system. The main peak position is consisten t 
with the cluster centroid estimated bv lPancino et all (2000). 

(ii) Both the RGB-a and the RGB-Mint populations have 
perturbed isodensity contour lines, showing structures sim- 
ilar to bubbles, shells and/or tails. While in the case of the 
RGB-a the complexity increases in the outer parts, where it 
is almost certainly due to low number statistics, in the case 
of the RGB-Mint we find complex structures in the central 
parts, where data points are more numerous. 

(iii) The RGB-Mint population is clearly elongated along 
the N-S direction in the inner parts, while in the external 
parts it seems to become elongated in the E-W direction. 
The main peak lies south of the RGB-MP peak, with a 
possible secondary peak north of it, that gives an evident 
asymmetric shape to the distribution. 

(iv) The RGB-a population is also elongated along the 
N-S direction in its inner parts; its main peak lies north of 
the RGB-MP peak. 

A simple monodimensional Kolmogorov-Smirnov test 
(K-S, see e.g.. iPress et alJll997l) as a function of a radial 
coordinate is not sufficient to properly assess the signifi- 
cance of these features. We thus used a two-dimensional 
generalization of the K-S statistical test, which was pro- 
posed and investigated with Monte-Carlo experiments by 
iFasano fe Franceschinil jl987l) , as a variant of an earlier idea 
b\ 1p^acoc1T(^83l) . This test, similarly to the usual K-S test, 
quantifies the probability P that two (two-dimentional) dis- 
tributions are extracted from the same parent distribution, 



using a more sophisticated definition of the maximum dif- 
ference D in the cumulative distributions. 

The results are summarized in Table The very low 
probabilities obtained in the comparison of the RGB-MP 
population with both the RGB-Mint and the RGB-a ensures 
us that their spatial distributions are significantly different, 
i.e., they cannot be drawn from the same parent distribu- 
tion. On the other hand, the probability obtained in the 
comparison between the RGB-Mint and the RGB-a popula- 
tions confirms that they are not significantly different from 
each other. 

We can also notice from FigureQthat the RGB-a pop- 
ulation appears more concentrated than the RGB-MP one, 
while we cannot say much about the RGB-Mint popula- 
tion, which has a complicated shape in its inner isodensity 
contour lines. A more quantitative evaluation confirms this 
impression: Table|2]reports the equivalent radii 2 of the three 
sub-populations, calculated with the ellipse parameters de- 
rived in Section 2] f° r the 0.9 and 0.6 isodensity contour 
lines. The RGB-a equivalent radii are ~2-3 times smaller 
than the RGB-MP correspondent radii. 



4 ELLIPTICITY 

To describe in a more quantitative way the shape and struc- 
ture of the sub-populations defined in Section we fitted 
ellipses to the isodensity contours of each sample. The spa- 
tial distributions were derived as in Section|3] with the same 
kernel, but using higher resolution grids of 20 and 50 pixel 
cells. We used the finer grid (20 pix) for the internal regions 
that are more densely populated, i.e., within the 0.6 isoden- 
sity level, and the coarser grid (50 pix) for the outer regions. 
As in Section [3] we defined densities in units of maximum 
(peak) density for each distribution, and we chose isodensity 
levels going from 90% to 10% of the peak density, in steps 
of 10%. Ellipses were fitted to each isodensity line, defined 
by those grid nodes that bracket the chosen density value 
(see Figure The best-fit ellipses are hereafter designated 
with the same notation used for the isodensity contours in 
Fig ure n i-e., 0.9 for the 90% ellipse and so on. 

Kho lopovl (119531) suggested that the best way to approx- 
imate equal density lines with ellipses in globular clusters is 



1 In what follows, it is useful to bear in mind that the WFI 
scale is 0.238 arcseconds per pixel. Thus, 100 pixels correspond 
approximately to 24 arcseconds. 



2 The equivalent radius for an ellipse of semi- axes a and b is 
defined as r e = \fab, i.e., the radius of a circle with the same 
area. 
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Figure 2. Example of the approximation of the 0.9 isodensity 
contour line with an ellipse, for the RGB-MP. Top Panel: The 
black dots represent the points in the 20 pixels grid that are closer 
to 90% of the peak value; the solid ellipse is the best fit. Bottom 
Panel Same as above, but this time with the 50 pixel grid: the 
final ellipse fit is slightly different. 



to use polar coordinates (r,ip). This method has the advan- 
tage that the deviations of points from the ellipse along the 
radial direction are close to the deviations in the direction 
perpendicular to the ellipse. In polar coordinates, the ellipses 
have the following form: 



= A + B ■ sin 2ij) + C ■ cos 2ip 



(1) 



The relations of the coefficients A, B and C with the 
usual ellipse parameters (Figure [3J are 
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Figure 3. Ellipse parameters notation: O is the ellipse centre; X 
and Y are the axes of the usual cartesian reference frame, while 
i/) and r are the angular and radial coordinates in the polar coor- 
dinate system; a and b are the ellipse semi-major and semi-minor 
axes, respectively; ip is the major axis inclination with respect to 
the X axis. 



As demonstra t ed by earlier studies 

l)Gever. Nelles fc Hood llQS.^ l. the "a priori" adoption 
of the cluster centre position can produce errors (i.e., an 
overestimation) on the ellipticity estimate. This point is 
even more important in our particular case, since we suspect 
that the centroids of the three sub-populations differ from 
each other. Therefore, we determined the coordinates of 
each ellipse centre as the mean coordinates of the input 
points on each isodensity line. The ellipse coefficients A, 
B and C where determined by least square approxim ation 
with singular value decomposition iPress et al 1 11997ft . An 
example of the result of the fitting procedure is shown in 
Figure H 



4.1 Ellipse centres 

Following the procedure described above, we determined the 
ellipse centres for each of the three sub-populations defined 
in Section [5] and for each isodensity level. The same proce- 
dure has been applied, for ease of comparison, also to the 
total RGB sample, defined as the union of the three sub- 
samples. The results are listed in TableEl the RGB-Mint and 
RGB-a centres are significantly different from the RGB-MP 
ce ntre. We can al s o com pare with the centre position found 
bv lPancino et all fcOOOft . which is (700±20,1900±20) in the 
WFI pixel system. As we can see, the centre of the global 
RGB sample is perfectly compatible with that estimate. 

A look at the centre position trend with the isoden- 
sity level (Figure shows that the RGB-MP, dominating 
the cluster population, has a quite stable centre position 
and coincide s reasonably we l l with the centre of ui Centauri 
measured bv lPancino et all (12000ft . except for its two inner 
isodensity levels. The RGB-Mint centre is slightly displaced 
to the east (~ 10") and significantly to the south (~ 1'), at 
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Table 3. The density peaks of the whole RGB sample and of the 
three sub-samples (in pixels), represented by the centres of the 
80% isodensity level fits. The distance of each population from 
the RGB total sample is calculated, both in pixels and in arcsec- 
onds. The last row reports the corresponding distance between 
the present RGB global samp le and the cluster centroid deter- 
mined bv lPancino et ahl l200CTl . 

Population Xc(pix) Yc(pix) a!(pix) d(") 



RGB-tot 


689±8 


1906±12 






RGB-MP 


701±9 


1954±17 


49±24 


12±6 


RGB-Mint 


651±10 


1671±16 


238±24 


57±6 


RGB-a 


761±9 


2194±8 


297±19 


71±5 


ui Centauri 


700±20 


1900±20 


12±32 


3±8 



least for the inner fits. This reflects the behaviour seen in 
Figure in the central parts the main peak of the popula- 
tion is clearly S-E of the RGB-MP population, while the be- 
haviour in the external parts becomes smoother. The RGB-a 
centre is instead slightly displaced to the west (~ 10"), and 
significantly to the north (~ 1'), at least for the few iso- 
density contour lines that we were able to fit with ellipses. 
Again, this reflects what seen in Figure We recall here 
that the literature estimates o f to Centauri's core radius go 
from a maximum of r c =2.89' llTra ger et a l.lll995l) to a min- 
imum of r c =1.4' llvan Leeuwen et aLB^OOoF ! so the observed 
displacements are comparable in size to the core radius of ui 
Centauri. 



4.2 Flattening and orientation 

As shown in Table the three sub-populations axial ra- 
tios have similar behaviours, within the errorbars. They do 
not show dramatic trends moving away from the centre, 
with both the RGB-MP and the RGB-Mint becoming slowly 
rounder away from their density peaks. For the RGB-a, due 
to the low number of objects in the external parts, only the 
inner isodensity levels could be fitted. The weighted aver- 
ages of the axial ratios shown in Table 0] are: < (b/a) >= 
0.81 ±0.01 for the RGB-MP, < (b/a) >= 0.81 ±0.06 for the 
RGB-Mint and < (b/a) >= 0.78 ± 0.11 for the RGB-a. 

To compare our results with previous work, we per- 
formed the same measurements on the entire RGB sample, 
resulting from the union of the three sub-samples defined 
in Section |21 In particular, Figure |5] shows the comparison 
with Geyer et al. (1983 - see their Table 4): the overall 
agreement is reasonably good, with a marginal discrepancy 
in the region between 2'. 5 and 4'. 5 from the cluster cen- 
tre. Moreover, our average axial ratio and ellipticity (e = 
1 - b/a) for the whole RGB sample, < (b/a) >= 0.89 ± 0.04 
and < e >= 0.11 ± 0.04, compare well with previous re - 
sults like < (b/a) > = 0.83 ± 0.03 iWhite fc ShawlHl987l) . 
< e >= 0.12 ± 0.04 jGever et al.lll"983r) . or < e >= 0.077 
llvan Leeuwen et alJl2000fl . 

Finally, while the flattenings of the three sub- 
populations appear rather similar, the orientations of the 
best-fit ellipses are instead quite different (see Table The 
RGB-MP major axis is always close to the E-W direction, 
with < <p >~ —4° ± 10, and its inclination seems to increase 
slowly with the distance from the cluster centre, although at 
these low flattenings the uncertainty on the inclination angle 
can be substantially higher than the formal errors quoted in 
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Figure 4. Displacement of the ellipse centres for the three-sub- 
populations. Ellipses are denominated along the abscissae in frac- 
tions of the peak value, as described in the text. The shaded area 
represents the cluster centre position determined by Pancino et 
al. (2000), (X c = 700 ± 20 pix, Y c = f900 ± 20 pix). The three 
sub-populations are marked with different symbols: filled dots for 
the RGB-MP, filled squares for the RGB-Mint and filled triangles 
for the RGB-a. The top panel shows the displacements in the E- 
W direction, i.e., along the X CCD axis, while the bottom panel 
shows the displacements along the N-S direction, i.e., the Y CCD 
axis. The errorbars are smaller than the symbols. 



Table B UGever et al.lll983T) . The RGB- Mint shows instead 
a complex structure in the central parts (Figure 0, with a 
possible double peak. But its 0.6-0.4 isodensity levels (the 
most reliable ones) are always inclined along the N-S direc- 
tion, with < ip >~ 100°, while the outer parts suddenly drop 
to < if >~ 25°, much closer to the E-W direction. Finally, 
the RGB-a is oriented along the N-S direction, although el- 
lipses could be fit only to the inner parts, with < tp > close 
to 90°-100°, except for the innermost fit. Again, these quan- 
titative estimates closely reflect what seen in Figure 



5 METALLICITY SEGREGATION 

Recentlv lJurcsikl jl998T) . using a compilation of good-quality 
spectroscopic data from the literature, has shown that the 
bright giants (V ^ 12.75) in u Centauri, belonging to the 
two metallicity groups with [Fe/H]^ —1.75 and [Fe/H]^ 
— 1.25, have a weird spatial distribution. Each of the two 
groups occupies one half of the cluster: no star of one group 
is found on the other side and the two centroids are sepa- 
rated by 6.2' (see her Figure 1). The separation line runs 
approximately perpendicular to the apparent minor axis of 
the cluster, which is roughly oriented towards the Galactic 
centre: the most metal rich giants are on the southern half, 
tha t faces the Galaxy. 

Ilkuta fc Arimota i2000f) questioned this effect. They 
used the same compilation of data and showed that no spa- 
tial segregation of stars with different metallicity is evident. 
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Table 4. Axial ratio and orientation for the best-fit ellipses for 
the whole RGB sample and for each sub-sample. Each best-fit 
ellipse is designated with the fraction of the peak density for each 
population, as described in the text. The axial ratio (b/a) and the 
inclination angle ip (i.e., the angle in degrees, counted from west 
to north) are shown. 




Population 


grid 


ellipse 


(b/a) ± 6(b/a) 


ip ± Sip 


IV Vj£3-TOb 


£AJ 


n Q 


n Q9 -i- n m 
u.yz m u.ui 


AO, 9-1-^7 

<-±y.z z\z o. I 




90 


U . o 


n qi -4- o ni 
u.yo nz u.ui 


9£ 1 4- A Q 




90 


n 7 

U. 1 


n Q1 -4- 09 
u.yo m u.uz 


1 7 X. -t- 1 X. 




20 


0.6 


0.93 ± 0.01 


19.5 ± 2.9 




50 


0.5 


0.93 ± 0.02 


18.0 ± 4.3 




50 


0.4 


0.93 ± 0.02 


20.9 ± 4.7 




50 


0.3 


0.92 ± 0.02 


31.9 ± 3.3 




50 


0.2 


0.89 ± 0.01 


32.0 ± 2.7 




50 


0.1 


0.84 ± 0.01 


3.7 ± 1.7 




50 


0.05 


0.84 ± 0.03 


6.0 ± 1.2 


DPR MP 

rx.Vjx £>- ivl ir 


90 
ZU 


Q 
U . I J 


H9 -1- 0^ 
U.oZ ZIL U.UO 


1 A -4- R 1 
-JLU."* It O.l 




20 


0.8 


0.80 ± 0.12 


3.8 ± 2.4 




20 


0.7 


0.79 ± 0.05 


7.0 ± 1.7 




20 


0.6 


0.81 ± 0.04 


7.0 ± 1.4 




50 


0.5 


0.84 ± 0.03 


12.4 ± 2.6 




50 


0.4 


0.85 ± 0.02 


15.8 ± 2.7 




50 


0.3 


0.88 ± 0.02 


22.7 ± 3.2 




50 


0.2 


0.90 ± 0.01 


41.9 ± 2.8 




50 


0.1 


0.85 ± 0.04 


8.6 ± 2.2 




50 


0.05 


0.85 ± 0.02 


10.9 ± 1.6 


RGB-Mint 


20 


0.9 


0.85 ± 0.16 


-8.0 ± 8.9 




20 


0.8 


0.73 ± 0.10 


-6.1 ± 2.4 




20 


0.6 


0.75 ± 0.02 


100.2 ± 0.8 




50 


0.5 


0.77 ± 0.03 


100.8 ± 1.6 




50 


0.4 


0.87 ± 0.03 


102.2 ± 2.8 




50 


0.3 


0.90 ± 0.02 


33.0 ± 6.4 




50 


0.2 


0.88 ± 0.02 


24.3 ± 3.8 




50 


0.1 


0.83 ± 0.06 


6.3 ± 2.2 


RGB-a 


20 


0.9 


0.66 ± 0.11 


47.2 ± 10.7 




20 


0.8 


0.86 ± 0.12 


97.9 ± 7.2 




20 


0.7 


0.81 ± 0.04 


112.7 ± 4.7 




20 


0.6 


0.83 ± 0.03 


108.7 ± 2.9 



Ho wever, we wou ld like to note here that, as clearly stated 
by IJurcsikl jl998fl . the effect is only visible if one applies 
the specified selections, i.e., only for bright stars of the two 
extreme met allicity groups. The effect was instead partially 
confirmed bv lHilker fc Richtlerl ll200(J) . who measured abun- 
dances of stars in to Centauri with Stromgren metallicity in- 
dexes. Their Figure 15 shows a clear segregation of the most 
metal-rich stars in the southern half of the cluster, but they 
could not confirm the segregation of the most metal-poor 

stars. 

■ ■ metal-rich group in IJurcsikl |l998), with [Fe/H]> 
— 1.25, is mainly a sub-sample of our RGB-Mint population, 
which has a pronounced peak just ~ 1' south of the cluster 
centre (see Figure and Table |^J, a value that is roughly 
compatible with Jurcsik estimate (~ 3'). We thus easily ex- 
plain the observed segregation of her metal-rich group, since 
when one draws randomly a sample of metal-rich stars, there 
is a higher probability to find them close to the main peak 
of the density distribution, i.e., in the southern half of the 
clus ter. 

IJurcsikl lll998l) metal-poor group, with [Fe/H]< -1.75, 
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Figure 5. Axial ratio for the whole RGB population (black dots), 
compared with the measures by Geyer et al. (1983), which is 
represented by a shaded area (see both their Table 4 and our 
Table EJ. 



is instead a sub-sample of our RGB-MP, or more precisely 
it represents the lowest metallicity tail of the RGB-MP. We 
tried to isolate this sub-population by dividing our RGB- 
MP sample vertically, in the (I,B — I) plane, in three sub- 
samples, equally wide in colour: (i) the bluest RGB-MP sub- 
sample, the RGB-MP1, most probably corresponding to Ju- 
rcsik's metal-poor group; (ii) the intermediate sub-sample, 
the RGB-MP2 and (Hi) the reddest sub-sample, the RGB- 
MP3. A two-dimensional KS test, like the one described in 
section gives a low probability that the RGB-MP 1 is ex- 
tracted from the same parent distribution of the RGB-MP2 
(P12 =1.42 1CT 6 ) or RGB-MP3 (Pi 3 =3.92 1(T 5 ). However, 
we were not able to measure any significant difference in 
the peaks positions or in the concentration of the RGB-MP 
sub-samples. 

We thus are unable to explain the me tal-poor part of the 
spatial metallicity segregation found by Ijjrrcjsfid l ) on 
the basis of our photometric catalogue. If lJurcsikl ( 199E ) ef- 
fect will be confirmed and understood in the future, it could 
mean that an additional, metal-poor sub-population exists 
in u Centauri, with its own distinct properties. Otherwise, 
the observed spatial segregation of that group of stars could 
simply be due to a sta tistical fluctuation, explaining why 
Hilkcr & Richtlcr (2000) were not able to confirm the metal- 
poor part of the segregation effect. More data on abundances 
of a significant sample of stars with [Fe/H]^ —1.75 are thus 
needed to completely explain this second half of the puzzle. 



6 SUMMARY AND DISCUSSION 

Using the nomenclature defined in lPancino et alJ i200fj) . we 
have selected three samples with B < 16: the RGB-MP, 
the RGB-Mint and the RGB-a. We have shown that the 
three samples are significantly different in their spatial and 
structural properties. In particular: 

(i) Both the RGB-Mint and the RGB-a main density 
peaks are shifted by ~1' with respect to the cluster centre. 

(ii) Both the RGB-a and the RGB-Mint populations are 
elongated on a direction that is perpendicular to the elon- 
gation of the RGB-MP population. 

(iii) The RGB-a population is significantly more concen- 
trated than the RGB-MP population. No firm conclusion 
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can be reached for the RGB-Mint, due to the complex iso- 
density contour line's shapes in its central parts. 

(iv) Both the RGB-Mint and the RGB-a show complex 
and perturbed isodensity contour lines, resembling bubbles, 
shells and tails. While for the outer isodensity contours this 
can be due to statistical undersampling (especially for the 
RGB-a), in the inner parts these features are probably real. 

(v) Finally, we were able to explain the metal-ric h part of 
the w eird spatial metall icity segregation, found by Ijurcsikl 
lll99Sft and confirmed bv lHilker fc Richtlerl (l2000l) . thanks to 
the peculiar surface density distribution of the RGB-Mint 
and, in particular, to its centre displacement with respect to 
the RGB-MP centre. 

Let us discuss these findings in the framework of the 
current scenarios for the formation and evolution of uj Cen- 
tauri. As briefly mentioned in Section 1, there is a consistent 
body of evidence supporting the fact that uj Centauri built 
up by itself at least part of the chemical elements that we 
can observe today. In a standard self-enrichment scenario, 
we would expect the more metal-rich (and younger) popula- 
tions to be more centrally concentrated than the metal-poor, 
dominant population. This is almost always the case for the 
dwarf galaxies of the local group, with only a few exceptions 
jHarbeck et alJl200ll) . and it is exactly what we observe for 
uj Centauri (point (iii) above). 

Unfortunately, while there is no doubt on the fact that 
self-enrichment must be one of the fundamental ingredients 
of a s uccessful formation and evolution theory fo r uj Cen- 
tauri jNorris fc Da Costal fl995l: ISmith et al.ll200Cl) . we are 
dealing with a complex set of observational facts, contain- 
ing conflicting evidence. In fact, a few of the observational 
properties of uj Centauri, concerning the structure, the shape 
and the kinematics (Section^ are not ea sily accommodated 
into a simple self-enrichment scenario Norris et alJ Il997t 
iPancino et all hoOO: Ferraro et al . 2002b), and suggest the 
possibility of a merger or accretion event. However, the sim- 
ple merging of two or more single metallicity clusters cannot 
account for the br o ad metallicity dist ribution of the RGB 
jNorris et alJll99rJ: Isinith et al.llboo'oTl and the high speed 
of ordinary, already formed globular clusters in the present 
potential well of t he Milky Way makes th is kind of merg- 
ing quite unlikely llThurl fc JohnstonlEo OZ*) 3 . The evidence 
presented here, concerning the structure and shape of the 
different sub-populations in the RGB of uj Centauri, sup- 
ports these findings and thus confirms the need for a more 
sophisticated scenario, that takes into account all the obser- 
vational evidence collected so far. 



6.1 A Complex Dynamical History 

Let us first discuss the relative orientation of the three RGB 
sub-populations. It is now well established that the elon- 
gated shape of the whole cluster (dominate d by the RGB- 
MP population) is mainly due to rotation l|Hardindll965t 
iMerritt fc Tremblavlll994i : iMevlan fc Mavorlll986l) . which is 
consistent with the picture of a dynamically young and not 
completely relaxed cluster: the cluster's peak rotational ve- 
locity is 7 km s _1 at 11 pc from the centre iMerritt et alJ 



3 In a lower potential environment, like in the Fornax dwarf 
spheroidal, the chance would be substantially higher. 



1997) and the relaxation times for uj Centauri are of the or- 
der of magnituide of its age, being l ogt re =9.73 yr in the 
core and \ogt r h =10 yr at half mass (lHarrislll99fi . 

Thus, it becomes tempting to explain the elongations 
of the RGB-Mint and RGB-a populations in terms of rota- 
tional velocity, too, but in this case it would be necessary to 
assume that these two populations rotate around a perpen- 
dicular axis with respect to the RGB-MP. An inspiring com- 
parison is posed by the recent work by Sarzi et al. (2000), 
who examine an example of galaxy that has undergone a 
major merging or accretion event in its past. The signature 
of such an event is the simultaneous presence of (1) an or- 
thogonally elongated bulge with respect to the disk, and of 
(2) two rotation curves, perpendicular to each other, one for 
the host galaxy and another for the accreted component. 

Now, we have shown here (point (ii) above) that the 
first signature could indeed be present in u) Centauri: we 
have (at least) two components with orthogonal elongations. 
What can we say about the rotation curve s? It has been 
demonstrated in the past jNorris et alJll997^) that while the 
metal-poor stars in uj Centauri rotate, no sign of rotation 
is evident for the metal-rich stars. Again, the metal-rich 
stars of Norris correspond mainly to our RGB-Mint pop- 
ulation, while too few data are presently available for the 
RGB-a population. So, the second signature is only par- 
tially present in uj Centauri. Clearly, to fully demonstrate 
or falsify the point, one needs more precise radial velocities 
for a much larger sample of RGB-Mint and RGB-a stars, 
and a model which is best suited to small systems embed- 
ded in a common potential well. For example, if we consider 
simulations o f the merger of two globular clusters with un- 
equal masses dMakino et al.llT99iri . the most massive object 
retains a larger shar e of th e initial orbital angular momen- 
tum, as lNorris et alJ dl997l) point out to explain why we see 
rotation for the metal-poor stars and not for the metal-rich 



6.2 An Accreted Component? 

Another interesting point has recent ly emerged from the 
coupling of the photometric cata log bvlPancino et al.l feOOCft 
and the proper motion work bv Ivan Leeuwen et al] J200fj) . 
i.e., that the RGB-a population appears to have a differ- 
ent bulk proper motion of | Sfj, | =0.8 mas yr -1 with respec t 
to the main population of uj Centauri llFerraro et al]l2002bl) . 
This evidence suggests that the RGB-a could be an accreted 
population, captured by the main body of uj Centauri. We 
thus should expect the RGB-a to have a different centre from 
the main cluster population, too, and it is exactly what we 
find here (point (i) above). Moreover, if a merging event has 
really taken place in the past history of uj Centauri, how long 
ago do we expect it to have happened? We already know that 
uj Centauri has a very long relaxation time, comparable to 
its age. We also know that the RGB-a is ~2-3 times more 
concentrated than the RGB-MP (point (iii) above), and the 
fact that it is almost one core radius away from the cen- 
tre suggests that it is moving in a significantly less dense 
environment, so it could have survived as a self-gravitating 
system for many Gyr. It is however puzzling that the (few) 
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radial velocities measured for the RGB-a so far are not so 
different from the cluster average 4 . 

So, a self-consistent set of observational facts, support- 
ing the occurrence of an accretion or merger event in the 
cluster past history, is beginning to take shape. We would 
like to note here that this fact in not necessarily in contradic- 
tion with the occurring of self-enrichment, nor with the hy- 
pothesis that ui Centauri is the remaining nucleus of a larger 
body (a dwarf galaxy), accreted and partially disrupted by 
the Milky Way. If we consider the case of the Sagittarius 
dwarf spheroidal (Sgr), which is the only one showing some 
resemblance with the case of u> Centauri, we see that more 
than one GC seems to be associated with the system. In par- 
ticula r, one of the se (Terzan 7) has a metallicity of [Fe/H]~- 
0.62 <HarriJl99d) compared with the [Fe/H]~-1.59 of M 54, 
which may be the nucleus of the Sgr. Interestingly, the RGB- 
MP has a similar metallicity to that of M 54, and the RGB-a 
to Terzan 7 5 . 

6.3 A Promising Framework 

An interesting scenario, that could accommodate all of the 
obse rvational eviden c e collected so far, was discussed by, 
e.g.. iFreemanl dl985l) . iNorris et all (Il997l) and ISmith et all 
(2000): the so-called merger within a fragment scena rio, de- 
scendi ng fro m the general framew ork proposed by ISearld 
<1977l) and ISearle fe Zrnnl lll97Sf) . In this framework, a 
conglomerate of star-forming sub-systems or regions could 
evolve within a large cloud and a common potential well 
(a fragment), each section evolving with slightly different 
timescales, and slightly different chemical properties. In this 
context, the chemical evolution of the RGB-MP and of the 
RGB-Mint could have been tightly related to each other, 
especially since there are reaso ns to believe that the RGB- 
Mlnt is younger by a few Gyr jHughes fc Wallersteinfe oOO ) 
and could have been enriched by the ejecta of the RGB-MP 
stars. 

Later, a dwarf galaxy with its own globular cluster sys- 
tem could form, and a few globular clusters (like the RGB-a, 
or even the RGB-Mint) could spiral towards the system cen- 
tre, while the remaining clusters could be stripped by the in- 
teraction with the Milky Way, alo ng with most of the dwar f 
galaxy halo. Recent calculations bv Bromm fe Clarke! ( 2002) 
supp ort this line of reasoning (see also iFellhauer fc Kroupal 
2002), together with the example of the Sagittarius dwarf 
galaxy discussed above. 

Although still speculative, this idea deserves to be fur- 
ther explored since it appears the most promising to explain 
ui Centauri's puzzling properties. 
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